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ABSTRACT 

We present interferometric measurements of the Sunyaev-Zel'dovich (SZ) effect towards the galaxy 
cluster Abell 370. These measurements, which directly probe the pressure of the cluster's gas, show the 
gas distribution to be strongly aspherical, as do the x-ray and gravitational lensing observations. We 
calculate the cluster's gas mass fraction in two ways. We first compare the gas mass derived from the SZ 
measurements to the lensing-derived gravitational mass near the critical lensing radius. We also calculate 
the gas mass fraction from the SZ data by deprojecting the three-dimensional gas density distribution 
and deriving the total mass under the assumption that the gas is in hydrostatic equilibrium (HSE). We 
test the assumptions in the HSE method by comparing the total cluster mass implied by the two methods 
and find that they agree within the errors of the measurement. We discuss the possible systematic errors 
in the gas mass fraction measurement and the constraints it places on the matter density parameter, 

Subject headings: cosmic background radiation — cosmology: observations, galaxies: 
clusters:individual(Abell 370)-techniques: interferometric 



1. INTRODUCTION 

Clusters of galaxies, by virtue of being the largest known 
virialized objects, are important probes of large scale 
structure and can be used to test cosmological models. 
Rich clusters are extremely massive, ~ 10 15 M Q , as indi- 
cated by the presence of strongly gravitationally lensed 
background galaxies and by the deep gravitational poten- 
tial necessary to explain both the large velocity dispersion 
(> 1000 km s~ x ) in the member galaxies and the high 
measured temperature (> 5 keV) of the ionized intraclus- 
ter gas. Dynamical mechanisms for segregating baryonic 
matter from dark matter on these mass scales are difficult 
to reconcile with observations and standard cosmological 
models, and so within the virial radius the mass compo- 
sition of clusters is expected to reflect the universal mass 
composition. Under the fair sample hypothesis, a cluster's 
gas mass fraction, which is a lower limit to the its baryonic 
mass fraction, is then a lower limit to the universal baryon 
fraction, i.e., f gas < f B . 

The luminous baryonic content of galaxy clusters is 
mainly contained in the gaseous intracluster medium 
(ICM). The gas mass is nearly an order of magnitude larger 
than the mass in optically observed galaxies (e.g., White 
et al. 1993, Forman & Jones 1982). Hence, the gas mass 
is not only a lower limit to the cluster's baryonic mass, it 
is a reasonable estimate of it. 

The intracluster medium has largely been studied 
through observations of its x-ray emission. The ICM is 
hot, with electron temperatures, T e , from ^5 to 15 keV; 
rarefied, with peak electron number densities of n e ~ 10 -3 



cm , and cools slowly (t coo i > t Hubble), mainly via ther- 
mal Brcmsstrahlung in the x-ray band. The x-ray sur- 
face brightness is proportional to the emission measure, 
S x oc J n^k(T e )dl, where the integration is along the line 
of sight, and so, under simplifying assumptions, the gas 
mass can be calculated from an x-ray image deprojection 
and the measured gas temperature. Since the sound cross- 
ing time of the cluster gas is much less than the dynamical 
time, one may reasonably assume that, in the absence of a 
recent merger, the cluster gas is relaxed in the cluster's po- 
tential. The total binding mass can be extracted from the 
gas density and temperature distribution under this as- 
sumption. A significant body of work exists in which the 
gas mass fraction, f g , is measured in this way, with f g mea- 
surements out to radii of 1 Mpc or more (White & Fabian 
1995; David et al. 1995; Neumann & Bohringer, 1997; 
Squires et al. 1997; Mohr, Mathiesen & Evrard 1999). 
The mean cluster gas mass fraction within approximately 
the virial radius was calculated in Mohr et al. (1999) to 
be (0.0749 ±0.0005)/i" 3 / 2 . Here, and throughout the pa- 
per, we assume the value of the Hubble constant to be 
Ho = 100ft. km s" 1 Mpc" 1 . 

To derive the gas mass fraction from x-ray imaging data, 
one is required to deproject the surface brightness into a 
model for the density distribution. As the x-ray emis- 
sion is proportional to the square of the gas density, the 
gas mass measurement can be biased by clumped, multi- 
phase gas, should it be present. Also, the emission from 
the cores of relaxed clusters may be dominated by cool- 
ing flows, which complicate the interpretation of the x-ray 
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data and may bias the result strongly if not taken into 
account (Allen 1998; Mohr et al. 1999). In addition, the 
x-ray surface brightness is diminished in proportion to its 
distance; S x oc 1/(1 + z) 4 , where z is the rcdshift of the 
cluster, neglecting experiment-specific K-corrections, and 
so it becomes increasingly difficult to make sensitive x-ray 
measurements of the ICM as the cluster redshift increases. 
We present a scheme for measuring the gas mass fractions 
with the Sunyaev-Zel'dovich effect which is different from 
the x-ray method in a number of ways, and also provides 
an independent measurement of f g . 

The Sunyaev-Zel'dovich (SZ) effect is a spectral distor- 
tion of Cosmic Microwave Background (CMB) radiation 
due to scattering of CMB photons by hot plasma (Sunyaev 
& Zel'dovich 1970). The SZ effect can be detected signifi- 
cantly in galaxy clusters, where the ionized intracluster gas 
serves as the scattering medium. A small fraction, < 1%, 
of CMB photons are inverse-Compton scattered and, on 
average, gain energy. At frequencies less than about 218 
GHz, the intensity of the CMB radiation is diminished 
as compared to the unscattered CMB, and the SZ effect is 
manifested as a brightness temperature decrement towards 
the cluster. This decrement, AT$z, has a magnitude pro- 
portional to the total number of scatterers, weighted by 
their associated temperature, ^ T ^ z ^ -f n eT e dl, where n e 
is the number density of electrons, T e is the electron tem- 
perature, Tqmb is the temperature of the CMB, and the 
integration is again along the line of sight. Note that this 
is simply proportional to the integrated electron pressure. 
Also, the magnitude of the SZ decrement is independent of 
rcdshift, so as the long as the cluster is resolved (the exper- 
iment's characteristic beamsize is not larger than the angle 
subtended by the cluster), the SZ effect can be measured 
towards arbitrarily distant clusters. 

A cluster's gas mass is directly proportional to its in- 
tegrated SZ effect if the gas is isothermal. So under 
the isothermality condition, an image deprojection is not 
strictly required to obtain the gas mass. The cluster's gas 
mass fraction can be calculated by comparing the inte- 
grated SZ decrement, in effect a surface gas mass, to the 
total cluster mass in the same volume. The total mass can 
be measured with strong or weak gravitational lensing, for 
example. The SZ images may also be dcprojcctcd to infer 
the three-dimensional gas mass and the HSE mass. Since 
the SZ decrement is directly propotional to the electron 
density, the SZ image deprojection will not be affected 
strongly by clumped gas. Thus, the cluster's gas mass 
fraction can then be measured as a function of cluster ra- 
dius as well. 

Recent cluster gas mass fraction measurements from SZ 
effect observations are presented in Myers et al. (1997). 
In this work, the integrated SZ effect is measured using 
a single radio dish operating at centimeter wavelengths. 
The integrated SZ effect is used to normalize a model for 
the gas density from published x-ray analyses, and this 
gas mass is compared to the published total masses to de- 
termine the gas mass fraction. For three nearby clusters, 
A2142, A2256 and the Coma cluster, Myers et al. find 
a gas mass fraction of (0.061 ± 0.011)/i _1 at radii of 1-1.5 
h^ 1 Mpc; for the cluster Abell 478, they report a gas mass 
fraction of (0.16 ± O.OW)^- 1 . 

In this work, we describe a method to calculate cluster 



gas mass fractions from intcrferometric SZ observations. 
Here, the shape parameters are derived directly from the 
SZ dataset rather than from an x-ray image. We apply this 
method to the SZ effect measurements towards the cluster 
Abell 370, which were made as part of an SZ survey of 
distant clusters, the first results of which were reported in 
Carlstrom, Joy & Grego (1996), and are further reported 
in Carlstrom et al. (1997). We choose this primarily be- 
cause it has been studied at optical and x-ray wavelengths, 
and so allows a comparison of the SZ data with other ob- 
servations. Upon detailed investigation, it is apparent that 
this cluster is one of the most difficult in our sample to an- 
alyze, with significant ellipticity and complicated optical 
and x-ray structure; as such, it serves as a test of the gas 
mass fraction analysis method, which we plan to use on a 
large sample of clusters. 

With the intcrferometric SZ measurements, comple- 
mented by observations at other wavelengths, we measure 
the cluster's gas mass fraction in two ways. First, we cal- 
culate the surface gas mass from the SZ measurements 
and the surface total mass from strong gravitational lens- 
ing observations and models. Second, we measure the gas 
mass fraction in a manner similar to x-ray analyses. The 
gas mass is inferred from a deprojected model; the total 
mass is determined from the spatial distribution of the 
gas under the isothermal hydrostatic equilibrium (HSE) 
assumption. We test the assumptions made in the HSE 
analysis by comparing the total cluster masses derived in 
the two methods. 

The optical and x-ray observations of this cluster are 
discussed in Section 2, and the SZ observations are dis- 
cussed in Section 3. The method for modeling the SZ 
data is presented in Section 4, and the gas mass fraction 
results and the systematic uncertainties are discussed in 
Section 5. The cosmological implications of the results 
and plans for future work are discussed in Section 6. 

2. OPTICAL AND X-RAY OBSERVATIONS 

The spatial distribution and redshifts of a number of 
Abell 370's constituent galaxies have been measured op- 
tically. Mellier et al. (1988), using their 29 best galaxy 
spectra, find Abell 370 has a mean redshift z = 0.374 
and velocity dispersion a v = 13401^° km s _1 . These 
observations also show that the cluster is dominated by 
two giant elliptical galaxies, at (aj2ooo = 02^ 39 m 52.7 s , 
<5j20oo = -01° 34' 20.3") and (a J20 oo = 02^ 39 m 53.2 s , 
<5j20oo = —01° 34' 57.6"), with a projected separation of 
about 40" in the north-south direction. The major axes 
of these two dominant galaxies are also oriented nearly 
north-south. The Mellier et al. optical isopleth map shows 
two density peaks, with similar north-south orientation, al- 
though these arc not centered on the two elliptical galaxies. 
The 29 galaxy velocities do not appear to be distributed in 
a Gaussian manner, although the statistics on this mod- 
est sample are not conclusive. Should this be significant, 
its may suggest that the cluster is bimodal or that it is 
not yet fully relaxed, and the assumption of a simple one- 
component model in hydrostatic equilibrium would not be 
appropriate. The mean velocities in the two isopleth peaks 
are not found to differ significantly, however, nor arc the 
calculated barycenters of the two peaks, and so it is diffi- 
cult to make a firm conclusion on the state of the cluster 
based on these data alone. 
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Abell 370 was the first cluster found to have an asso- 
ciated strongly gravitationally lensed arc (z ar c — 0.725; 
Soucail et al. 1988). In addition to this giant arc, which 
apparently consists of multiple images of a single source, 
three faint pairs of images, each pair presumed to be lensed 
images of a single galaxy (Kneib et al. 1993), and a radial 
arc (Smail et al. 1996) have been found to be associated 
with this cluster. Gravitational lens models have been 
used to map the total mass distribution in the core region 
of Abell 370. Kneib et al. (1993) find the position and 
elongation of the giant arc and the positions of the multi- 
ple images are best reproduced by a bimodal distribution 
of matter, with the two centers of mass nearly coincident 
with the cluster's two large elliptical galaxies. A subse- 
quently discovered radial arc is consistent with this model 
without modifications (Smail et al. 1996). 

Abell 370 was observed with the High Resolution Imager 
(HRI) of the ROSATx-my satellite. In the fiOS^T band- 
pass (0.1-2.4 kcV), x-ray Brcmsstrahlung emission from 
a plasma at gas temperatures typical of clusters depends 
only weakly on the gas temperature. In Figure 1, we show 
an HRI image of A370, the result of <~ 30 kiloseconds 
of observation time. The bright source to the north, at 
"72000 - 02 h 39 m 55 s , <5 J20 oo = -01° 32' 33", is coin- 
cident with what appears in the Digital Sky Survey to 
be a nearby elliptical galaxy, NPM1G -01.0096 (Klemola, 
Jones & Hanson 1987), and does not appear to be asso- 
ciated with A370. The image indicates that the intra- 
cluster medium also is elongated in the north-south direc- 
tion, and the emission is strongly peaked at the position 
of the southern dominant elliptical galaxy, numbered 35 
in Mellier et al. (1988), and is offset from the center of the 
large scale emission. The cluster gas distribution does not 
appear clearly bimodal at this resolution. 

We determine the x-ray surface brightness from the 
ROSAT HRI image, which has been filtered to include 
only PHA channels 1-7 in order to reduce the background 
level (David et al. 1997) and blocked into 8"x 8" pixels. 
A circular region with a radius of 24" centered on the 
pointlikc x-ray source to the north of the cluster was ex- 
cluded from our modeling. We fit a /3-model to the data; 
this model described in more detail in Section 4, which 
parametrizes the x-ray surface brightness with a core ra- 
dius, 9 C , and a power law index, (3. The surface brightness 
of a spheroidally symmetric, isothermal gas distribution 
will have elliptical contours, and so we also fit for the axis 
ratio. 

We find the best-fit axis ratio of the x-ray image to 
be 0.73±0.03, with major axis exactly north-south, and 
the center position to be aj2ooo = 02^ 39 m 53 s , <5j2ooo = 
—01° 34' 33". The best fitting shape parameters for the 
elliptical model are (3= 0.70±0.03 and 6> c = 85.0±3.5". If 
the axis ratio is set to equal 1.0, the best fitting shape 
parameters are (3= 0.72±0.03 , 6 C = (70.1 ± 3.2)" (90% 
confidence on a single parameter). 

Abell 370 was also observed with the Einstein IPC for 
4000 seconds. In the IPC image, the northern point source 
is not apparent. This may be because the FWHM of the 
instrument point spread function is about 90" at the mid- 
dle of the 0.4-4.0 keV spectral range of the instrument. 
(It improves to about 60" at higher energies.) The point 
source is offset by less than 120" from the center of the 



extended emission, and so may not be resolved in this im- 
age. It may also be explained by source variation; if the 
point source is an AGN, its surface brightness may have 
significantly varied between the 1979 IPC observation and 
the 1994-1995 HRI observations. The image looks less el- 
liptical than docs the HRI image, but because the point 
source emission cannot be excluded from the cluster emis- 
sion in the IPC data, it is not useful for constraining the 
spatial distribution of the cluster gas. 

Abell 370 has also been observed with the ASCA x- 
ray satellite. Using 37.5 kiloseconds of ASCA GIS/SIS 
observations, Mushotzky & Scharf (1997) determined the 
emission-weighted average temperature of the cluster to 
be 7.13+ 

1x83 keV. Analyzing the same data, Ota, Mitsuda, 
& Fukazawa (1998) find the emission-weighted tempera- 
ture to be 6.6to 9 keV. Both of these values are system- 
atically lower than the value obtained in the first analysis 
of these data by Bautz et al. (1994), who obtain a value 
of 8.8 ± 0.8 keV. The discrepancy is presumably because 
the Mushotzky et al. and Ota et al. analyses used the val- 
ues of the response matrices of the ASCA telescope and 
x-ray detectors which have been refined since the Bautz et 
al. analysis (M. Bautz, private communication). For our 
analysis, we adopt the Ota et al. value of 6.6t|xg keV, be- 
cause this published analysis includes detailed discussion 
of the spectral fitting procedure. 

The bright point source evident in the R OS A T HRI im- 
age has not been removed from the ASCA spectra in any 
of the three published analyses. If this emission originates 
in Bremsstrahlung emission from hot gas in the ellipti- 
cal galaxy, its characteristic temperature is expected to be 
less than 1 or 2 keV (Sarazin & White 1988 and references 
therein). Since the sensitivity of ASCA is optimized for 
higher energies than this, we would not expect the mea- 
sured temperature to be greatly contaminated. We assess 
whether the source is extended, and find that it is pointlike 
at the resolution of the HRI with a 2.4 a certainty. As the 
source is unresolved, this emission is also consistent with 
a galactic cooling flow; if this were the case, the tempera- 
ture again shouldn't be greatly contaminated. Should the 
emission originate from an AGN in the galaxy, though, 
the effect on the measured temperature is more difficult 
to predict, and will depend strongly on the index of the 
power law of the AGN spectrum. 

Without further measurements, we cannot rule out that 
the cluster's spectrum is contaminated by emission from 
an AGN. It is also possible that if the cluster is in 
fact bimodal with two distinct subclusters, the measured 
emission-weighted average T e will be a value between the 
temperatures of the two subclusters. Lacking more defini- 
tive measurements, we continue with our assumption that 
the gas is isothermal at kT — 6.6lg'g keV. The estimated 
effects of the described uncertainties on the results are dis- 
cussed in Section 5.3. 

We compare A370's measured line-of-sight velocity dis- 
persion with our adopted gas temperature, and find 
/3 spec = \xm v a^jkT e = 1.87lJ'5 9 . ^ the energy per unit 
mass contained in the galaxies is equal to that in the gas, 
and the cluster is spherical, we expect (3 spec to be ~ 1. If 
the temperature measurement is biased low, e.g., from an 
AGN, this will contribute to the /3 spec discrepancy. Elon- 
gation of the cluster potential along the line of sight or a 
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alignment of two (or more) subclusters will also enhance 
the measured f3 spec . 

The optical and x-ray observations suggest that A370 is 
not a relaxed cluster, and may be heavily substructured, 
but due to insufficient signal-to-noise ratios and contam- 
inating sources, these observations still leave some ambi- 
guity. 

3. SZ OBSERVATIONS 

We observed Abell 370 for 50 hours over nine days with 
the Owens Valley Radio Observatory (OVRO) Millimeter 
Array in 1996 August and 43 hours over six days with the 
Bcrkeley-Illinois-Maryland Association (BIMA) Millime- 
ter Observatory in 1997 June-August. We set our pointing 
center to a J20 oo = 02 71 39 m 52.5 s , <5 J20 oo = -01° 34' 20", 
the position of the northern dominant elliptical galaxy. 

We outfitted the millimeter telescopes with centimeter- 
wavelength receivers to broaden the instrument's angu- 
lar resolution to the large angular scales typical of clus- 
ters. The receivers are based on low noise High Electron 
Mobility Transistor (HEMT) amplifiers (see Pospieszalski 
1995), operating from 26 to 36 GHz. The single-mode re- 
ceivers are constructed to respond only to circularly polar- 
ized light. The receivers interface with each array's delay 
lines and correlator. The standard observational scheme 
entailed interleaving 20 minute cluster observations with 
observations of a strong reference source near the clus- 
ter, which allows monitoring of the instrumental amplitude 
and phase gain. The amplitude gain drifts were minimal; 
the gain changes were less than a percent over a many hour 
track, and the average gains were quite stable from day 
to day. The amplitude gains were calibrated by compar- 
ing the measured flux of Mars with that predicted by the 
Rudy model for Mars' whole disk brightness temperature 
(Rudy, 1987) at the observed frequency and the apparent 
size of Mars as indicated by the Astronomical Almanac. 
The Rudy model is a radiative transfer model with an es- 
timated accuracy of 4% (90% confidence) at centimeter 
wavelengths. Data taken when the projected baseline was 
within 3% of the shadowing limit were not used. 

The OVRO Millimeter array consists of six 10.4 me- 
ter telescopes. The OVRO continuum correlator was used 
to correlate two 1 GHz bands centered at 28.5 GHz and 
30 GHz. The system temperatures ranged from 40-60 K, 
depending on elevation and atmospheric water content. 
The primary beams of the telescopes were measured holo- 
graphically and can be approximated as Gaussian with a 
full width at half maximum (FWHM) of 252". The data 
were calibrated and edited using the MMA data reduction 
package (Scoville et al. 1993) taking care to remove data 
taken during poor weather or which show any anomalous 
phase jumps. 

We image the data using the DIFMAP package (Shep- 
ard, Pearson, & Taylor 1994). Examining images of A370 
made with projected baselines greater than 2.0 kX, we find 
a point source 45" to the cast of the map center (map cen- 
ter is the pointing center) in both continuum bands. Its 
measured flux density in the 28.5 GHz band, attenuated by 
the instrumental primary beam pattern, is +0.69 ± 0.10 
mJy and is +0.84 ± 0.10 mJy in the 30 GHz band. This 
source, 0237-0147, is discussed further in Cooray et al. 
1998. We combine the two continuum channels, and find 
the best fitting model is a point source with flux density of 



0.725 + 0.07 mJy. The point source model is removed from 
the u-v or spatial frequency, data set in order to construct 
the SZ decrement map. A Gaussian taper with FWHM of 
1.2 k A is applied in the u-v plane to the combined dataset 
and map was CLEANed, restricting the CLEAN compo- 
nents to the central 200", about the size of the decrement 
in the primary. The map presented in Figure 2a is made 
with a restoring beam of Gaussian FWHM 59" x 86", with 
contour intervals of 1.5 a. The rms noise in the map is 50 
/z Jy beam -1 , or 14.9 /zK, and the integrated flux density 
of the source is —1.54 ± 0.17 mJy. 

We employed nine of the 6.1 meter telescopes of 
the BIMA array for SZ observations, using the same 
centimeter- wave receivers we used at OVRO. We center 
the 800 MHz bandwidth of the BIMA digital correlator 
at 28.5 GHz. We achieved system temperatures of 30-55 
K, scaled to above the atmosphere, depending on eleva- 
tion and atmospheric water content. Holographic antenna 
pattern measurements were also made with this system. 
The primary beam at BIMA is nearly Gaussian with a 
396" FWHM. The data were edited and reduced using 
the MIRIAD data reduction program (Sault, Teuben & 
Wright 1995), taking care to the remove any spurious in- 
terference from the spectral channels and dropping the low 
signal-to-noise channels at the spectral filter edges. 

Again using DIFMAP, we find the same point source 
at about 45" to the east of the pointing center in a high 
resolution (projected baselines longer than 1.4 kX) image. 
The point source measured at BIMA in 1997 has flux den- 
sity of +0.70 ± 0.17 mJy. The primary beam attenuation 
at the point source position in the BIMA system is ~ 6% 
less than the attenuation at this position at OVRO, and 
so the point source observations are consistent with the 
flux being constant in time. We model and subtract the 
point source from the data set, apply a Gaussian taper 
with FWHM of 0.8 fcA in the u-v plane and construct a 
map of the decrement. The map was CLEANed, restrict- 
ing the CLEAN components to the central ~ 300" of the 
image, about the size of the primary beam. The restoring 
beam used to make the map in Figure 2b. is a Gaussian 
with FWHM 91" x 95", and the contours presented are at 
intervals of 1.5 a. The rms noise in the map is 180 /iJy 
beam -1 , or 31.4 ^K, and the integrated flux density of the 
source is —4.34 ± 0.52 mJy. 

It is instructive to compare the general characteristics of 
the SZ images to the x-ray map. At the resolution of both 
the OVRO and BIMA instruments, the gas is extended in 
the north-south direction, though less markedly than in 
the x-ray map. The substructure suggested by the lensing 
and velocity dispersion data is not evident here. We pro- 
duce maps from each of the datasets, increasing the image 
resolution, but the gas still does not show significant sub- 
structure or bimodality. The SZ decrement distribution is 
not peaked at the pointing center, but south of it, at a 
position about halfway between the two elliptical galaxies 
(this will be quantified in the next section). The detailed 
appearance of the SZ map, especially the shapes of the 
least significant contours, depends on the specific method 
of point source removal and the CLEANing of the data, 
and so it is not useful to compare the spatially filtered and 
CLEANed SZ and x-ray images at such a level of detail. 
As we will discuss in the following section, the quantitative 
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analysis of the SZ data is done using the u-v data directly. 
4. MODELING 

In order to assess the data more quantitatively, we fit 
a parametrized model to the SZ brightness distribution. 
Since the SZ data are taken in the spatial frequency, or 
u-v, domain, we do the model comparisons in u-v domain 
as well. 

We measure the SZ temperature decrement in units of 
antenna temperature, T a , which relates an observed in- 
tensity change to a Rayleigh- Jeans temperature change, 
AJ„ = ^-AT a . In order to recover the true blackbody 
temperature decrement, the ATjjj measured and reported 
here should be multiplied by a factor of 1.021, as we are 
not strictly in the Rayleigh- Jeans limit. The SZ temper- 
ature decrement, Tcmb , is proportional to the Compton 
y-parameter, 



y : 



kcfx 



n e {l)T e {l)dl, 



(1) 



where k is Boltzmann's constant, <j T is the Thomson scat- 
tering cross section, m e is the electron mass, n e is the elec- 
tron density, T e is the electron gas temperature, and the 
integral extends along the line of sight (dl). The propor- 
tionality depends on the observing frequency; it depends 
also on the electron temperature when relativistic correc- 
tions are included 6 (Rephaeli 1995, Challinor & Lasenby 
1998). 



At 28.5 GHz, 



TcMB 



T.92 y in the non- relativistic 



Rayleigh- Jeans approximation, where we adopt the value 
of the CMB of Fixsen et al. (1996) derived from the COBE 
FIRAS measurements, Tqmb = 2.728 K. Including the rel- 
ativistic corrections for kT e — 6.6 keV, ^ Tsz = —1.86 v. 

' T C MB y 

We fit a simple model to the SZ data. The /3-model 
(Cavaliere & Fusco-Femiano 1976, 1978) is frequently used 
to fit the density profiles of galaxy clusters. The spheri- 
cally symmetric, isothermal /3-model describes n e , the den- 
sity distribution of the gas, as varying as a function of 
cluster radius, r: 

/ r 2\ -3/3/2 

n e (r) = n eo ( 1 + -» . (2) 



where n eo is the number density of the gas at the center of 
the spheroidally symmetric gas distribution; r c , the core 
radius, is a characteristic size of the cluster; and /3 is the 
power law index. 

Since the cluster appears elliptical in projection, we gen- 
eralize the spherically-symmetric /3-model and allow the 
gas density profile to be spheroidally symmetric, i.e., with 
biaxial symmetry. In the spheroidally symmetric model, 
the electron number density in a prolate spheroid is a func- 
tion of rj 2 = r 2 + a 2 z 2 , where r and z arc the radial and 
height coordinates in the cylindrical coordinate system, 
and a < 1 is the axis ratio of the two unique axes. If 
the spheroid is oblate, then the density is a function of 
rj 2 = a 2 r 2 + z 2 , with a < 1. The electron density then 
follows the distribution: 



1 + r 2 



-3/3/2 



(3) 



If the cluster gas is isothermal and its symmetry axis is in 
the plane of the sky, this electron density distribution leads 
to the following two-dimensional SZ temperature decre- 
ment: 

_M 



AT (0) = AT(0) 




(4) 



where 9 V = tj/Da, Da is the angular diameter distance, 
8 C = r c /DA 7 and AT(0) is the temperature decrement at 
zero projected radius, 

AT(0) oc T e n e0 [ ( 1 + ( -) \ dl, (5) 



where the integral, dl is along the line of sight. Formally, 
this integral extends from the observer along the line of 
sight through the cluster infinitely; in practice, a cutoff 
radius for the cluster is used. The /3-model distribution 
of the electron density projects to a /3-model distribution 
of the SZ decrement if the system obeys spherical or ellip- 
soidal symmetry. 

We perform a x 2 analysis of the /3-model, by compar- 
ing /^-models to the combined BIMA and OVRO datasets. 
We vary the following parameters: centroid position, /3, 
9 C , axis ratio (defined to be < 1), position angle (defined 
counter clockwise from north), and AT(0). The position 
and flux density of the radio-bright point source are also 
fit. The fitting procedure is conducted in several steps. 
The model corresponding to each set of the seven clus- 
ter fit parameters and the three point source parameters 
is multiplied by the primary beam response. The Fourier 
transform of the result is compared directly with the in- 
terferometer data. We use the holographically determined 
primary beams when modeling the data, and the entire 
datasets are used to do the analysis. The inner u-v ra- 
dius cutoff is determined by the shadowing limit, the limit 
where one telescope would partially block another, i.e., 
when the projected baseline is less than the diameter of a 
telescope dish. For the BIMA data this limit is 0.58 kX 
and for the OVRO data it is 1 kX. The \ 2 statistic is 
minimized, and the best fit values are determined using a 
downhill simplex method. 

Performing the fitting procedure, we find the best fit 
parameter values for A370 are: 9 C = 93", (3 = 1.77, AT(0) 
= —609 [iK , and axis ratio = 0.64 with major axis exactly 
north-south. The best fit central position is 20" to the 
south of the pointing center, at aj20oo = 02^ 39 m 53.2 s , 
<5j20oo = —01° 34' 40.4", about halfway between the two 
giant elliptical galaxies, and 7" north of the fitted centroid 
for the x-ray image. The x 2 statistic for the best fit values 
is 135798 for 136136 degrees of freedom, yielding a reduced 
X 2 of 0.9975. As points of comparison, we fit to the data 
a null model for the cluster (with the point source compo- 
nent), which gives a reduced x 2 of 0.9986; we also fit the 
cluster to a negative point source, finding a reduced x 2 of 
0.9986. 

In order to determine the uncertainties in these fitted 
parameters, the x 2 statistic is derived for a large range 



"The change in spectral intensity due to the Sunyacv-Zcl'dovich effect is calculated in Equation (4-8) Challinor & Lasenby (1998): 



A(T SZ /T) RJ 
f(x) = 3.58. 



w[a:coth(a;/2) — 4 + 8 e f(x)], where x 



hv 



and i 



The last term corrects for relativistic effects. At 28.5 GHz, 
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of [3, 9 Cl and AT(0), keeping centroid position, position 
angle, and point source position and flux fixed to the best 
fit values. Presented in Figure 3a is a contour plot of the 
(3 and 9 C fit results, where AT(0) is allowed to assume its 
best fit value at every pair of (3 and 9 C . The axis ratio is 
left to assume its best fit value at each set of (3, C , and 
AT(0); it varies from 0.60 to 0.68 for points with A* 2 < 
5 from the best fit point. The centroid position, when left 
to assume its best fit value at each point, does not ap- 
preciably change, varying less than 5". Figure 3b shows 
the contours derived when fitting the data to a spherical 
model for the gas distribution, i.e., the axis ratio is fixed 
to a value of 1.0; in this case, the best fit parameter values 
are 6 C = 34", (3 = 0.86, AT(0) = -785 fj,K, the reduced 
X 2 statistic was 0.9976. 

The full line contours are marked for A% 2 = 2.3, 4.61, 
and 6.17 which indicate 68.3%, 90.0%, and 95.4% confi- 
dence, respectively, for the two-parameter fit. The dashed 
lines indicate Ax 2 = 1.0, 2.71, and 6.63; the projection 
onto the (3 or 8 C axis of the interval contained by these 
contours indicate the 68.3%, 90% and 99% confidence in- 
terval on the single parameter. The values of AT(0) for 
which A% 2 < 1 at each (f3, 9 C ) point range about 15% from 
the best fit value. This illustrates that (3 and 9 C are corre- 
lated strongly and arc not individually well constrained 
by these data. This poor constraint is at least partly 
due to the low declination of A370, which makes sampling 
non-redundant spatial frequencies difficult at OVRO and 
BIMA. The shape parameters fit from the x-ray data are 
not consistent with those from the SZ data. This could be 
a result of shocks and complexity in the gas phase which 
affect the x-ray emission and SZ effect differently. This 
discrepancy is difficult to resolve with the information at 
hand; for the purpose of the next section, investigating the 
gas mass fraction of the cluster using the SZ effect, we use 
the SZ fitted parameters only. 

We also fit the data with a /3-model which includes a 
truncation of the gas distribution at a given radius. With 
truncation radii from 300" to 1000", the shape parameters 
(/?, C , AT(0)) best fit values do not change appreciably, 
much less than the variation contained within the A\ 2 < 
1 region. The \ 2 statistic changes only minimally for dif- 
ferent cutoff radii, with A% 2 less than 0.1. This indicates 
that the possible systematic uncertainty in the fitted pa- 
rameters due to modeling the data without a cutoff is not 
significant. 

Since the optical and x-ray data suggest that Abell 370 
may have a bimodal gravitational potential, we also fit the 
data with a pair of circular /3-models, each allowed inde- 
pendent shape parameters and position. The best fit two- 
component model has one component centered at ~ 50" 
south of the pointing center, 10" south of the large south- 
ern elliptical galaxy; the best fit position of the second 
component is ~ 5" north of the pointing center, i.e., 5" 
north of the large northern elliptical. The southern com- 
ponent has best-fit parameters (3 = 0.65, 9 C = 10", and 
AT(0) = —765 /iK; the northern component has best-fit 
parameters (3 = 0.83, 9 C = 27", and AT(0) = -605 /iK. 
Adding a second component reduces the \ 2 statistic by 2 
when 3 new parameters are introduced (compared to the 
single (3- model with axis ratio and position angle free), and 
therefore is not a significantly better fit to the data. Al- 



lowing the axis ratios of the two components to vary does 
not improve the fit. 

5. ANALYSIS 

5.1. Comparison of SZ Gas Mass to the Strong Lensing 
Mass Estimate 

As indicated by Equation 1, the SZ brightness at any 
point in the two-dimensional projected image is simply 
proportional to the integrated electron density along the 
line of sight, if the gas is isothermal. Under the isothermal 
assumption, we can directly measure the total number of 
electrons in the gas contained in the cylindrical volume of 
a chosen radius, the long axis of which is defined by the 
line of sight. The total mass in the ionized phase can be 
calculated from this assuming a value for the number of nu- 
cleons per electron. If the gas has solar metallicity, as mea- 
sured by Anders & Grevesse (1989), the nucleon/electron 
ratio is 1.16. The nucleon/electron ratio changes less than 
1% for values of the metallicity from 0.1 to 1.0. 

The interferometric measurements recover much of the 
total SZ decrement on the angular scales measured, i.e. 
the integrated flux from the data comprises 40-50% (de- 
pending on the instrument) of the integrated flux in a best 
fit model to these data. And so the unmodeled interferom- 
eter data provides a strong lower limit to the integrated 
SZ decrement at the angular scales of interest. The model 
is fitted to the u-v data in order to estimate the full, two- 
dimensional decrement, and thence the surface gas mass. 
This method does not assume anything about the state of 
the ICM other than that it is isothermal. 

We fit the SZ data in the three-dimensional parameter 
space 0.4 < (3 < 4.0, 10" < r c <250", -2000 /iK < AT(0) 
< —200 /iK, both for elliptical and circular models. At 
each (3, 8 C , and AT(0) point, we calculate the surface gas 
mass and the \ 2 statistic, from which confidence limits for 
the surface gas mass are determined. Although the 68% 
confidence region contains a large range of (3 and 9 C val- 
ues within 68% confidence, as was evident in Figure 3, the 
gas mass, which depends on all three gridded parameters, 
is constrained relatively well. In Figure 4, we show the 
derived surface gas mass at radius 65" as it varies with (3 
and 9 C in the circular /3-model fitting. The gas mass in 
this geometry is 1.51^6 x 10 13 ft~ 2 M Q . The temperature 
decrement is allowed to assume its best fit value at each 
point for this figure, although for the quantitative analy- 
sis, the full range of AT(0) is used. The isomass surfaces, 
shown in greyscale, follow the shape of the confidence in- 
terval contours, indicating that sets of parameters which 
fit the data well will predict the same gas mass. This is to 
be expected, since under the isothermal assumption, the 
SZ flux is directly proportional to the gas mass, and the 
fit parameters must reproduce the same observed flux at 
the angular scales where the flux is best measured. 

We derive a surface gas mass to compare with the lens- 
ing mass using our fitted model parameters. The lenscd 
arc in Abell 370 has a radius of curvature of about 30" 
centered nearly halfway between the dominant galaxies; 
this center is at about the same position of the gas den- 
sity centroid in the one component [3 model. We calcu- 
late the surface gas mass in a cylindrical volume centered 
at this position and with an elliptical (axis ratio = 0.64) 
cross-section, using the fits to the elliptical /3-model. We 
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choose a major radius of 40" since the lens model should 
be most accurate near the lensed arc's radius. The SZ- 
derived surface gas mass with the elliptical cross-section is 
5.4±}; 2 x lQ 12 h- 2 M Ql at 68% confidence. Included in the 
error estimates are the uncertainty due to the fit parame- 
ters, the SZ absolute calibration uncertainty, and the gas 
temperature measurement uncertainty. The uncertainty is 
not appreciably smaller if we restrict the shape parame- 
ters to those which typically found in x-ray image analyses 
of clusters, i.e., (3 < 1.5. The uncertainties of the masses 
derived from the elliptical fits are smaller than those from 
the spherical fits because the confidence contours follow 
the isomass contours more closely. 

We also derive the gas mass in the model with two com- 
ponents. The best fitting two-component /3-model param- 
eters are integrated in a cylinder with radius 40", cen- 
tered on the centroid of the best fit single component 
model. This yields a gas mass of 4.84 x 10 12 /i~ 2 M Q , consis- 
tent with the single component model. Since the bimodal 
model doubles the parameter space over which fits must be 
made, making a comprehensive parameter fit unfeasable, 
we approximate the statistical uncertainty from the fit to 
be the same as in the single model fit, about 20%. 

These surface gas masses are compared to the cluster's 
total mass in the same volume implied by the strong gravi- 
tational lensing measurements. The surface total mass can 
be inferred from a model of the cluster mass distribution 
which predicts the observed gravitational lensing. We cal- 
culate the gravitational mass using the lensing model in 
Kncib et al. (1993) for the same volume, centered at the 
SZ model fit center. Although the uncertainties on the 
model parameters are less than ten percent, we make a 
conservative estimate of 20% for the uncertainty, to al- 
low for variations of this model. Comparing this mass, 
9.61^9 x l0 13 h~ 1 M Ql to the single-component gas mass 
yields a gas mass fraction, f g , of (0.0561q'q^5)/i _1 ; the two- 
component model yields a slightly lower gas mass fraction 
of (0.046 ± 0.013)/!- 1 . The cluster's angular diameter dis- 
tance was calculated assuming Hm = 0.3 and fl\ = 0. If 
Om = 1-0, £!a = 0, the angular diameter distance to A370 
will be 6.5% smaller, as will f g . If Q,m — 0.3 and £!a = 
0.7, f g will be - 9% larger. 

5.2. Comparison of SZ Gas Mass to Hydrostatic, 
Isothermal Mass Estimates 

Since strong gravitational lensing in clusters is relatively 
rare and is restricted to the cores of clusters, and weak 
lensing analyses are published for only selected clusters, 
we also consider the more general means of calculating the 
cluster gas fraction, deprojection of the density model and 
the HSE assumption. We assume the gas is in hydrostatic 
equilibrium, is isothermal, and is spheroidally symmetric; 
for simplicity, we assume the symmetry axis is in the plane 
of the sky. The surfaces of constant electron density are 
then concentric ellipsoids. We then compare the gas mass 
and the HSE mass from the deprojected model to deter- 
mine the gas mass fraction. 

Specifically, to determine the gas mass, we extract the 
central electron density, n eo , from the deprojected /3-model 
and measured electron temperature by performing the in- 
tegral along the line of sight in Equation 5 and integrating 
Equation 2. We do this for each set of (3, 9 C , and AT(0). 



These calculations require an assumption about the geom- 
etry of the cluster, since the core radius and extent of the 
cluster along the line of sight arc not known. For an oblate 
ellipsoid, the axis of symmetry is the cluster's minor axis, 
and the core radius in the line-of-sight direction is equal to 
the cluster's observed major axis; for a prolate ellipsoid, 
the axis of symmetry is the major axis, and the core radius 
in the line-of-sight direction is equal to the minor axis. 

The fit parameters are then used to constrain the total 
mass. Hydrostatic equilibrium implies 



1 



■Vp. 



gas 



(6) 



where p gas and p gas are the gas density and pressure, re- 
spectively. The cluster's gravitational potential, $, can 
be related to the total mass density, p gra v, by Poisson's 
equation, 



V 2 $ = 4nGp 



grav ; 



(7) 



where G is the gravitational constant. We can solve for 
the density of the cluster's gravitational mass by combin- 
ing Equations 6 and 7: 



Pgrav 



AttG 



(8) 



We relate the pressure, density, and temperature of the 
gas through the equation of state: 



Pgas — 



Pgas^Te 
flTTlp 



(9) 



where k is Boltzmann's constant, \i is the mean molecu- 
lar weight of the gas, and m p is the proton's mass. To 
calculate /i, we again assume the gas has the solar metal- 
licity of Anders and Grevesse (1989) and that p is constant 
throughout the gas. Making the assumption that the gas 
is isothermal, we write Equation 8 in the form: 



Pgrav 



A-KGprrip 



V 2 lnp 9 



(10) 



Note that the gravitational mass density depends only on 
the shape of the gas distribution, and so is independent 
of the value of the central gas density and the gas mass 
fraction. Using the derived shape parameters, (3, 9 C , and 
projected axis ratio, a, and the measured gas temperature, 
we estimate the total mass density of the cluster. We again 
choose the simplest geometries, that of oblate and prolate 
ellipsoids, and integrate the density within the same vol- 
ume as we do the gas density. 

For the same range of f3, 8 C , and AT(0) as we used in the 
cylindrical geometry analysis, we calculate the cluster's el- 
lipsoidal gas mass, HSE mass, and gas mass fraction for 
both prolate and oblate geometries, and the A% 2 from the 
best fit parameters. From A% 2 = 1 range of fit param- 
eters, we derive the 68% confidence intervals for the gas 
mass and the mean gas mass fraction calculated within 
different major axis radii. 

We prefer to measure the masses and mass fractions 
in the largest volume permitted by our method, since 
the fair sample assumption is best at large radii and the 
cores of clusters may be affected significantly by physical 
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processes not included in our HSE model (cooling flows, 
galaxy winds, magnetic fields) . The largest scales on which 
we make our calculations are determined by the shortest 
baselines on which we detect the SZ effect. We calculate 
the statistical uncertainties in the f g measurement due to 
the shape parameter uncertainties on a number of scales, 
from 10" to 150". There is a broad minimum in uncer- 
tainty around radius 65". 

We calculate the gas mass fraction for the oblate and 
prolate spheroids at semi-major axis 65" (~ 330/i _1 kpc). 
The gas mass is the same for both oblate and prolate ge- 
ometries; the change in central density and volume when 
the geometry changes exactly compensate. For the oblate 
ellipsoid, the gas mass fraction is (0.064t^o24)^ _1 ; f° r the 
prolate ellipsoid, the gas mass fraction is (O.O^^Q^^h" 1 . 
We also calculate f g in a spherical volume of radius 65", 
using the fits and uncertainties from the spherical fits and 
find f g = (0.0801^4! )h~ 1 . The calculated values are com- 
piled in Table 1. 

We also calculate the gas mass fraction for both com- 
ponents of the bimodal model, using a gas temperature of 
6.6 keV for each and find that the gas mass fraction for 
each component is <~ 0.040/i _1 . 

The dark matter needed to produce concentrically el- 
lipsoidal isopotential surfaces has a significantly more as- 
pherical distribution than the potential. For the observed 
axis ratio of 0.64, this dark matter distribution is unphys- 
ical beyond a few core radii in the direction of the unique 
axis, as it requires the dark matter density to be nega- 
tive. This also implies that the gas mass fraction will vary 
spatially. The difficulty of this model supporting very el- 
liptical gas density distributions may be suggesting that 
the cluster is bimodal, but it is more likely an indictment 
of the simplified models we have used. We have considered 
the class of potentials with isopotential surfaces following 
concentric ellipsoids because they adequately describe the 
data and because their deprojections are straightforward. 
Such simplified models are seriously deficient, however, for 
the hydrostatic analysis. To produce a cluster with an ob- 
served axis ratio of less than 1 / y/2 in this formalism, the 
cluster mass must be partially comprised of dark matter 
with negative density. It is improper, then, to measure the 
dark matter with this method out past a few core radii for 
a highly elliptical cluster. 

The effect of using this simple ellipsoidal model to cal- 
culate f g is assessed in the following way. We construct a 
model cluster which has a physically motivated mass struc- 
ture, predict its observed SZ effect, and then attempt to 
recover the simulated cluster's mass using the same obser- 
vation and analysis protocol used for the true observations. 
We arrange the dark matter of the simulated cluster in con- 
centrically spheroidal shells with a /3-model profile. The 
dark matter's axis ratio is set to 2.0:3.5. Isothermal gas at 
7 keV is added in hydrostatic equilibrium with the cluster 
potential, and a simulated two-dimensional SZ decrement 
map constructed, assuming the cluster is at redshift z = 
0.37. The decrement map is sampled with the u-v coverage 
of a typical interfcrometric observation. Noise typical of 
a <~ 40 hour observation is added. The resulting u-v data 
are fit in the same manner described in Section 4, with 
the data fit to a concentrically ellipsoidal /3-model, and 
the gas mass and gas mass fraction determined with the 



HSE method. Two simulations are used, one model clus- 
ter is an oblate spheroid and the other a prolate spheroid, 
both with the symmetry axis in the plane of the sky. In 
both cases, the best fit axis ratio of the gas distribution 
was 0.79. These gas masses and gas mass fractions are 
compared to the actual model values. These are shown in 
Table 2. 

The model appears to be adequate at small radii, but 
deteriorates noticeably by a radius of 100". This suggests 
it is reasonable to measure the gas mass fraction with this 
model near the same radius we constrain the data well. 
We restrict our analysis of the real observations to within 
radius of 65"(^ 330/i~ 1 kpc), a region within which the 
approximation is valid. 

5.3. Systematic Uncertainties 

5.3.1. Comparing the Lensing and HSE Masses 

If the assumptions made in the hydrostatic isothermal 
analysis are valid, the HSE predicted mass should equal 
the cluster's lensing mass. We integrate the total mass 
density in Equation 10 in the same cylindrical volume in 
which the lensing mass is calculated. This comparison is 
potentially a means to discriminate between oblate and 
prolate models for the cluster, but, as discussed in Sec- 
tion 5.2, our model for ellipsoidal clusters breaks down at 
large radii. For the spherical model, this mass is calcu- 
lated for each set of shape parameters and the 68% confi- 
dence limits are derived. Using an electron temperature of 
6.6jZQ 9 keV, the spherical geometry HSE mass integrated 
in the 40" radius cylindrical volume is lA7to.H x 10 14 M Q , 
consistent with the lensing mass in the same geometry of 
1.59±3;|g x 10 li h- 1 M Q . The two-component /3-model pre- 
dicts about 1.25 x 10 14 M Q , and so it is also a consistent 
model. 

This agreement does not guarantee that the assumptions 
in the HSE analysis are valid, however. For this reason, we 
examine possible systematic effects. Some comparisons of 
cluster total masses derived by the HSE method to lens- 
ing masses suggest that the HSE method systematically 
underprcdicts the cluster's mass (e.g., Miralda-Escude & 
Babul 1994, Loeb & Mao 1994, Wu & Fang 1997). Some 
of the explanations suggested for the discrepancy include 
cluster ellipticity, non-thermal pressure support of the gas 
in the cluster core, multi-phase gas, and temperature gra- 
dients in the gas. In an examination of a large sample of 
clusters, however, Allen (1998) suggests such discrepancies 
can often be resolved by taking account of cooling flows 
when analyzing x-ray data, and ensuring that the lensing 
and x-ray masses probe the same line of sight. The mass 
discrepancy does remain for clusters in the Allen sample 
with small cooling flows or none, perhaps suggesting these 
clusters have undergone recent dynamical activity which 
have disrupted any pre-existing cooling flow. Such activity 
might invalidate the HSE assumption. 

We have been careful to ensure that the strong lens- 
ing and SZ models probe the same lines of sight. There 
is no evidence for a cooling flow in Abell 370, although 
there may be other concerns about using the measured 
emission-weighted gas temperature. 

5.3.2. Contamination of Emission- Weighted Temperature 
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If the measured average temperature is in error, i.e., due 
to contamination from a nearby AGN, the mass measure- 
ments will also be in error. The SZ gas mass is inversely 
proportional to the assumed temperature and the HSE 
mass is directly proportional to the temperature. The gas 
fraction from the HSE method is quite sensitive to tem- 
perature, fg OC 1/Tg. 

5.3.3. Polytropic Temperature Gradient 

An unresolved temperature gradient in the gas may sys- 
tematically affect the gas and HSE masses. If such a gra- 
dient is present, the true temperature in the central region 
may be higher than the emission- weighted temperature we 
use, and the fitted shape parameters from the isothermal 
SZ analysis may no longer accurately describe the density 
distribution. 

However, if the temperature of the intracluster medium 
declines slowly, and does not change appreciably over the 
angular scales to which we are sensitive, the interferomet- 
ric measurement of the gas mass fraction at these angular 
scales will in fact not be strongly affected. 

Currently, there are no strong observational constraints 
on temperature structure in moderately distant and dis- 
tant clusters, as there have been no suitable telescope fa- 
cilities for the task. However, an attempt to quantify tem- 
perature structure observed in nearby clusters is presented 
in Markevitch (1996). In this work, a slow decline in tem- 
perature with radius is observed, with the temperature 
falling to half its central value at 6-10 core radii. This 
structure may be approximately described by a gas with a 
polytropic index of 7 = 1.2. (For discussion of polytropic 
indices, see Sarazin (1988) and the references within.) If 
there temperature structure in Abell 370 and it is of the 
moderate variety presented in Markevitch (1996), it will 
not be a strong source of systematic uncertainty. 

Abell 370 has a long observation scheduled with the 
Chandra observatory. Chandra has the necessary spatial 
and spectral resolution to remove the effects of contam- 
inating sources in the field and with a long observation, 
could measure the spatial variation of the cluster temper- 
ature, should it exist. 

5.3.4. Inclination Angle 

The assumption that the cluster is biaxially symmetric 
with symmetry axis in the plane of the sky is certainly 
a simplification. Here, we estimate the effect on the gas 
mass fraction of an inclined symmetry axis. Analytical 
relationships between the inclination angle and the appar- 
ent shape of an biaxial ellipsoid are derived in Fabricant, 
Rybicki, & Gorenstein (1984). We derive the cluster gas 
mass and HSE mass for a cluster which reproduces the 
observed SZ map, but has an inclined symmetry axis. A 
biaxial cluster, when inclined, will retain a /3-model distri- 
bution with the same value of /?, but the central density 
and intrinsic axis ratio will change. The inclination angle 
i is measured between the symmetry axis and the line of 
sight; i = 90° for the analysis of Section 5.2. 

As inclination angle increases, so does the intrinsic axis 
ratio. 7 In order to preserve the observed SZ effect, The 

7 The relationship between intrinsic axis ratio, a^, and the inclina 
the observed axis ratio. For a prolate ellipsoid, = a— „ sln o 7 



value of the central density must vary inversely with the 
axis ratio. The gas mass, then, remains the same for all 
inclination angles. This is expected, as the gas mass is 
proportional to the total SZ flux, independent of its spa- 
tial arrangement. We evaluate the gas and HSE masses as 
a function of i, using the best fit parameters from A370. 
The gas mass fraction changes, but it changes significantly 
over a relatively small range of allowed inclination angles 
(see Figure 5). 

6. DISCUSSION 

6.1. Gas Mass Fractions at r 50 o 

To compare the f g we have measured within a fixed 
angular radius to f g measurements in clusters with differ- 
ent sizes and redshift, we extrapolate our measured f g to 
a fiducial radius. The ICM in nearby clusters has been 
observed to be distributed more uniformly than the dark 
matter (e.g., David et al. 1995), which is to be expected 
if energy has been added to the intracluster medium be- 
fore collapse or from galactic winds. If this is generally 
true, the gas mass fraction measured depends on the ra- 
dius within which the measurement is made. As suggested 
in Evrard (1997) and Metzler, Evrard, & Navarro (1998), 
we choose this radius to be that within which the aver- 
age density of the cluster is 500 times the critical den- 
sity, p c — H^y- These numerical simulations suggest 
that within this radius, t^qq, the cluster's baryon fraction 
should closely reflect the universal baryon fraction, if the 
current physical models of hierarchical structure formation 
are correct. We use the analytical expression of Evrard 
(1997) which describes the expected variation of f g with 
overdensity. This variation is found to be consistent with 
the f g variation reported in the David et al. (1995) sample. 

UrMTe)) = f 9 (rx)(^^)\ (11) 

where r\ — 0.17, f g (r 500 (T e )) is the gas mass fraction at 
rsoo, and rx is the radius within which the gas mass frac- 
tion is measured. We modify Evrard's expression for r^oo, 
derived for low redshift clusters, to include the change in 
the value of p c with redshift; p c {z) — p c (z — 0)(H/H o ) 2 , 
where H 2 = H 2 [(1 + zfn M + (1 + z) 2 (l - Q M - SI a) + Ov] 

f rp X 1/2 

, 500 (T e ) = (1.24 ± 0.09) ( 1QkeV(W ) /^Mpc. 

(12) 

The gas mass fraction values at rsoo, estimated from 
those measured at 65", are summarized in Table 1. 

This experiment best measures f g at a given angular 
scale, which corresponds in Abell 370 to an overdensity of 
~ 5000p c . This is not the optimal radius at which to com- 
pare with numerical simulations, since resolution is limited 
in the cores of the clusters, and the gas in the core may 
also be sensitive to additional physics not yet included in 
the models, e.g., magnetic fields and cooling. For these 
reasons, the corrections should be taken with some cau- 
tion. 

angle, i, for an oblate ellipsoid is a; = a ^ 1-a — si ^/'^ ^ 1 where a is 
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6.2. Constraints on Qm from f g 

Under the fair sample hypothesis, A370's gas mass frac- 
tion within rsoo, a lower limit to the cluster baryon frac- 
tion, should reflect the universal baryon fraction: 



fg < SB = 



(13) 



where f B is the cluster's baryon fraction, SIm is the ra- 
tio of the total mass density to the critical mass density, 
and SIb is the ratio of baryon mass density in the universe 
to the critical mass density. The cluster gas mass frac- 
tion measurements can then be used within the Big Bang 
Nucleosynthesis (BBN) paradigm to constrain Qm- 



Qm < ^b/ fg 



(14) 



The value of is constrained by BBN calculations 
and the measurements of light element abundances. The 
relative abundance of deuterium and hydrogen provides 
a particularly strong constraint on the baryonic matter 
density. A firm upper limit to Qb is set by the presence 
of deuterium in the local interstellar medium. This con- 
strains the value of fis to be less than 0.031/i~ 2 (Linsky et 
al. 1995). Measurements of the D/H ratio in metal-poor 
Lyman-a absorption line systems in high-rcdshift quasars 
put a tighter constraint on the baryonic mass density. Such 
measurements made by Buries & Tytler (1998) predict a 
value of n B = (0.019 ± 0.002)/i- 2 at 95% confidence. 

The gas mass fractions measured for A370 from both 
lensing and HSE methods range from 5 — 13%/i~ 1 . We 
consider the simplest f g measurement, that in the spher- 
ical model, and compare this gas mass fraction at r§oo 
to the Buries & Tytler (1998) value for Cl B . This gives 
an upper limit to the matter density parameter, Om < 
0.19^q' 1 o' 1_1 ) a t 68% confidence. However, the bimodal 
model gives a surface gas mass fraction at angular radius 
40" of 0.048/i _1 , a value which permits VLm to be as high 
as 0.40/i _1 in this scheme. (We note again that the de- 
pendence of f g through the angular diameter distance is 
weak at this redshift, with a change in f g of 5-10% when 
a wide range of cosmological parameters is used.) These 
values are consistent with the limits on Qm from obser- 
vations of supernovae, which are derived from geometrical 
arguments, rather than the Qb/ fg ratio. Depending on 
the method used to calibrate the sample, for a spatially 
flat universe, Garnavich et al. (1998) find Qm < 0.4 — 0.5 
at 68% confidence. 

6.3. Conclusion & Future Work 

We have measured the Sunyaev-Zel'dovich effect in the 
galaxy cluster Abell 370 and present spatially filtered im- 
ages from these data. The optical and x-ray observations 
of this cluster show a complicated and perhaps bimodal 
mass distribution. The SZ effect image, however, looks 
smoothly distributed and significantly aspherical. We have 
fit both one-component and two-component /3-modcls to 
the data and find that the two-component model does not 
fit significantly better. For both models, we calculate the 
gas mass fraction for the cluster using measurements of the 
total cluster mass from both the gravitational lens model 
(the "surface" gas mass fraction) and from the hydrostatic 
equilibrium assumption. When integrated in the same vol- 
ume, the HSE masses are consistent with the mass derived 



from the gravitational lensing model for both the one- and 
two-component models, lending support to the HSE as- 
sumption. The surface gas mass fraction measurement is 
made within an angular radius of 40" and the HSE gas 
mass fraction is made within a radius of 65". The gas mass 
fraction near the virial radius is derived from the gas mass 
fractions at 65" using a correction factor derived from nu- 
merical simulations. For the range of methods and models 
used, we find gas mass fraction values of <~ (5 — 13)h~ 1 %. 

Constraints on the Hubble parameter H can in princi- 
ple be derived from the SZ and x-ray measurements of a 
cluster. The SZ and x-ray observables depend on differ- 
ent moments of the electron density, and so the charac- 
teristic length scale of the cluster along the line of sight 
can be measured and the angular diameter distance of 
the cluster inferred. This is useful not only as a distance 
ladder-independent measurement of H , but, when com- 
pared with other H a measurements, can be used to explore 
possible systematic effects in the f g calculation, e.g., to 
constrain the deprojection of the gas distribution. How- 
ever, the quality of the x-ray imaging data in this case and 
the apparent disagreement between the SZ and x-ray fit- 
ted models do not permit putting a strong constraint on 
the Hubble constant. A long observation with the Chan- 
dra x-ray observatory towards this cluster is planned, and 
should help resolve these issues. 

The value of the baryonic mass fraction in any one clus- 
ter will be susceptible to systematic uncertainties which 
may be difficult to estimate. In order to use cluster gas 
mass fractions as a cosmological tool, one wants to ame- 
liorate the effect of these errors by studying a large sample 
of clusters. A SZ effect survey in galaxy clusters is being 
carried out by this group at the BIMA and OVRO obser- 
vatories, and an analysis of the gas mass fraction of the 
sample is in preparation. 
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Table 1 

Gas Mass Fractions for Abell 370 

geometry f ± S -£g\ f g (S 5ao ) 

cylinder (radius=40") 
cylinder (radius=40"), bimodal 
oblate ellipsoid (radius=65") 
prolate ellipsoid (radius=65") 
sphere (radius=65") 



(o.oselg^)*- 


-i 






(o-Mfitg;™)*- 


-i 






(0.06418:81)*- 


-i 


1.22 




(o.io6i8:8S)^- 


-i 


1.22 


(0.12918:8s)*- 1 


(o.o8o±g:8**)/»- 


-i 


1.22 


(0.09818-.8S)*- 1 



Note. — The surface f g (cylindrical geometry) is derived from the elliptical (3- 
model (with a = 0.64) and the two-component model for the SZ effect and the 
total mass from the strong gravitational lensing model of Kncib ct al. (1993). The 
ellipsoidal gas mass fractions are calculated from deprojections of both the elliptical 
and spherical /3-model fits to the SZ data and the isothermal HSE assumption. The 
major axis indicated in the first column and symmetry axes are assumed to be in 
the plane of the sky. The radius within which f g is measured is indicated in the first 
column; f g within r 50 o is estimated in Column 4, using a correction factor (Column 
3) from Equation 11. 



Table 2 

Comparison of Derived to True Gas Mass and Gas Mass Fraction 

Within 65" Within 100" 

Oblate Cluster, i = 90° 

derived value/true value derived value/true value 

total mass 1-05±^| 1.42±g;^| 

gas mass 0.98±^:^| 0.96±^2 
gas mass fraction 

Prolate Cluster, i = 90° 

derived value/true value derived value/true value 

total mass Q-SS-Vil 1.16±oil 

gas mass 0.8718:8? 0.88l8:?4 

gas mass fraction l-H-o^s °- 79 ^8:20 



0.08 u - au -0.12 

gas mass fraction 0.95l8:i8 °- 65 -8'.06 



Note. — We test the effect of using an ellipsoidally symmetric /3-model for the 
gas density distribution. We simulate SZ observations of an oblate cluster and a 
prolate cluster and analyze them with the methods of Sections 4 and 5.2. We 
compare the quantities we derive to the true simulated cluster's values. 
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A370: ROSAT HRI X-ray Image 




-01°36'00" ' 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 — I 

2 h 40 m s 57 s 54 s 51 s 48 s 

RA (J2000) 

Fig. 1. — X-ray image of the galaxy cluster Abell 370 observed with the ROSAT HRI. The x-ray image was filtered to include only PHA 
channels 1-7, binned in 8" x 8" pixels, and smoothed with an 10" Gaussian filter. The background level is 2.8 x 10 — 3 HRI counts sec -1 
arcmin ~ 2 and the countour levels are at 5.8, 7.2, 8.7, 9.5, 10.1, and 11.6 x 10 — 3 HRI counts sec -1 arcmin ~ 2 . The x-ray emission peak is 
associated with a pointlike source located ~ 2.5' to the north of the cluster center. 
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Fig. 2. — a. Image of Abcll 370 from OVRO observations. The restoring beam, shown in the lower left hand corner, has a Gaussian FWHM 
of 59" x 86". The RMS noise in the map is 50 fijy beam -1 , or 14.9 fiK, and the integrated flux density of the source is —1.54 ± 0.17 mjy. 
b. Image of A370 from BIMA observations. The restoring beam has FWHM of 91" x 95". The RMS noise in the map is 180 [ijy beam -1 , 
or 31.4 /iK, and the integrated flux density of the source is —4.34 ± 0.52 mjy. 
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Fig. 3. — Confidence regions from the joint fit to the BIMA and OVRO SZ data using a. an elliptical /9-modcl and b. a spherical /3-model. 
In both plots, the best fit (ft, C ) point is indicated with an x and the full line contours show the A% 2 = 2.3, 4.61, and 6.17 regions, which 
indicate 68.3%, 90.0%, and 95.4% confidence, respectively, for the two-parameter fit. The dashed lines show A\ 2 =1.0, 2.71, and 6.63 regions; 
the projection of these regions onto the (i or C axis indicate the 68.3%, 90% and 99% confidence interval on the single parameter. The 
centroid and point source flux and position are fixed at the best fit values, and AT(0) is allowed to assume its best fit value at each (/?, C ) 
point. 
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Gas Fraction vs. Inclination Angle 




Fig. 5. — The ellipsoidal gas mass fraction, calculated with the isothermal HSE method, as a function of inclination angle, i. The mass 
calculated using the best fit parameters from the elliptical /3-modcl fit. 



